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EXECUTIVE  SUMMARY 


This  report  describes  our  procedure  and  initial  study  of  the  feasibility  and  practicality  of  using  a 
frequency  modulated  (FM)  chirp  to  perform  a  fine  positional  localization  of  hydrophone  elements  of 
a  passive  array.  Results  are  analyzed  with  respect  to  variation  based  on  direction  of  the  sound  source 
relative  to  the  array  and  the  ratio  of  sample  rate  to  array  frequency.  Results  show  that  the  measure¬ 
ments  are  not  consistent  enough  that  a  single  measurement  can  be  taken  and  used  with  great  confi¬ 
dence.  However,  averages  of  about  half  a  dozen  measurements  appear  to  have  sufficient  consistency. 

Extraneous  sources  of  acoustic  energy,  which  were  not  planned  or  controlled  in  the  test  setup,  are 
shown  to  be  of  comparable  utility  in  localization  of  the  array  elements. 
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INTRODUCTION 


The  objective  of  this  document  is  to  address  the  feasibility  and  practicality  of  using  an  FM  chirp 
to  determine,  on  a  fine  scale,  the  location  of  individual  elements  of  a  linear  underwater  acoustic 
array. 

The  principle  was  proven  to  be  effective  in  controlled  conditions  based  on  a  test  conducted  in  Jan¬ 
uary  1994  .  In  this  earlier  test,  an  eight-element  array  was  used  with  similar  signals  in  approxi¬ 
mately  40  feet  of  water  in  San  Diego  Bay.  This  document  relates  to  the  first  attempt  to  use  this 
method  in  water  that  is  deeper  than  diver  depth  (so  the  straightness  of  the  array  cannot  be  directly 
assured  or  verified).  Since  signals  spanning  a  longer  time  were  used  in  this  test,  the  results  are  not 
directly  comparable  with  the  earlier  test.  However,  the  signals  used  in  both  tests  had  the  same  band¬ 
width  and  should  achieve  the  same  resolution. 

TEST  SETUP  RELEVANT  HARDWARE 

The  two  arrays  used  in  this  test  consisted  of  a  compound  array  (nesting  a  640-Hz  array  within  a 
160-Hz  array)  and  a  simple  160-Hz  array  separated  by  1  km  center  to  center.  The  total  number  of 
elements  in  the  two  arrays  was  originally  to  have  been  64. 

Digitization  and  telemetry  of  the  data  from  each  hydrophone  was  accomplished  in  eight  nodes, 
each  handling  data  from  eight  hydrophones.  Acoustic  data  were  digitized  to  14-bits  at  a  sample  rate 
of  3276.8  Hz.  The  first  hydrophone  from  each  node  was  sampled  simultaneously,  the  second  hydro¬ 
phone  from  each  node  was  sampled  12  psec  later,  and  each  node’s  third  hydrophone  was  sampled 
12  psec  after  that,  etc.  See  figure  1  for  a  sampling  timing  diagram.  In  this  way,  the  timing  skew 
between  channels  (hydrophones)  is  known,  and  has  a  maximum  value  of  7  X  12  psec  =  84  psec.  In 
terms  of  the  sample  period,  l  =  1/3276.8  =  305  psec,  so  the  maximum  skew  is  about  27.5%  of  a  sam¬ 
ple  period. 

Prior  to  the  deployment,  six  of  the  hydrophones  were  eliminated  due  to  manufacturing  problems, 
and  the  spacing  was  adjusted  accordingly.  Upon  deployment,  certain  other  hydrophones  stopped 
working  or  existed  with  reduced  functionality.  For  this  reason,  there  are  a  few  ways  in  which  the 
hydrophones  can  reasonably  be  numbered.  Refer  to  figure  2  for  an  explicit  diagram  of  the  arrays  and 
interelement  spacing.  The  primary  means  of  labeling  that  shall  be  used  is  sequential  numbering  of 
the  58  physically  present  hydrophones,  starting  at  the  shore  end. 

The  arrays  were  laid  along  the  bottom  as  straight  as  possible  along  a  line  of  bearing  of  135°  true, 
to  the  southwest  of  Point  Loma,  San  Diego,  CA.  The  water  depth  in  this  area  is  approximately 
90  meters.  See  figure  3. 

Two  active  sound  sources  were  also  deployed  on  the  ocean  floor.  Initially,  one  was  placed  normal 
to  the  high-frequency  array,  and  one  was  placed  in-line  between  the  arrays.  The  one  normal  was 
later  placed  at  a  location  normal  to  the  low-frequency  array.  The  fine  layout  of  the  arrays  and 
sources  is  shown  in  figure  4.  These  sources  emitted  the  FM  signals  at  scheduled  times.  Two  differ¬ 
ent  chirps  were  used.  The  high-bandwidth  chirp  was  0.25-second  long,  and  frequency  dropped  lin¬ 
early  from  1280  Hz  to  280  Hz.  The  low-bandwidth  chirp  was  also  0.25-second  long,  and  frequency 
dropped  linearly  from  640  Hz  to  140  Hz.  The  high-bandwidth  chirp  was  transmitted  at  the  top  of  the 
minute,  and  the  low-bandwidth  at  40  seconds  after.  Timing  accuracy  was  achieved  by  a  calibration 
with  a  Global  Positioning  System  (GPS)  receiver  before  the  source  was  deployed. 
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Timing  for  each  of  eight  telemetry  nodes.  Phone  1  from  each  node  is  simultaneous. 

Figure  1.  Timing  diagram. 
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sequential  numbering  of  15  phones  in  LF  array 
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Figure  4.  Approximate  source  locations. 


The  recording  system  was  implemented  in  a  486  computer.  It  was  capable  of  recording  up  to 
2.5  seconds  of  data,  beginning  at  a  particular  time  specified  as  HH:MM:SS  coordinated  universal 
time.  The  recording  computer  had  a  full-time  link  to  a  GPS  receiver.  All  data  used  for  this  report 
were  recorded  within  a  2-hour  period  on  the  afternoon  of  14  September  1994. 

PROCEDURE 

The  relative  delays  between  hydrophones  were  estimated  using  a  seven-step  process  that  is 
detailed  in  this  section.  An  eighth  step  is  also  shown,  which  converts  delays  to  angle  of  arrival. 

1,  Pare  down  the  matrix. 

The  data  was  initially  stored  as  a  matrix  A  of  N  X  64,  where  N  is  the  number  of  samples,  usually 
between  6000  and  8000.  Since  the  sound  sources  were  always  within  about  1000  meters  of  the 
arrays,  only  a  few  thousand  samples  were  necessary  to  acquire  the  signals,  and  the  later  samples 
were  eliminated.  Also,  the  data  matrix  was  64  columns  wide,  corresponding  to  the  original  64 
hydrophones  on  the  string.  The  21st,  28th,  29th,  37th,  44th,  and  60th  columns  were  eliminated, 
since  these  correspond  to  hydrophones  that  were  no  longer  physically  present  when  the  array  was 
deployed.  This  leaves  the  raw  data  in  a  N  X  58  matrix,  with  N  varying  on  a  file  by  file  basis  from 
1536  to  4096.  This  matrix  was  B  such  that, 

B  =  Ay  j^{21,28,29,37,44,60} 

2,  Zero  out  the  bad  channels. 

As  seen  in  figure  NO  TAG,  there  are  seven  hydrophones  that  were  not  producing  any  useful  sig¬ 
nal.  So  they  would  not  influence  the  statistics,  these  channels  were  replaced  with  zeros.  At  the  same 
time  phone  #40  (of  58)  was  inverted,  which  compensates  for  its  being  wired  180°  out  of  phase. 

The  channels  were  zeroed  instead  of  being  removed  because  they  are  useful  to  retain  as  place¬ 
holders.  An  array  with  uniform  spacing  can  still  be  formed  with  the  inclusion  of  these  virtual 
phones. 

These  operations  can  be  easily  accomplished  by  an  element-by-element  multiplication  of  row  of 
Bby  a  “correction  vector”  c. 

C  =  [0  1  1  1  1  1  1  0  1  1...1  011110001111111101  -1  1...1] 

Bj  =  Bj  •  c  ,  Vi 

3,  Replica  correlation. 

The  intended  method  of  postprocessing  was  to  filter  the  received  signal  using  a  matched  filter 
based  on  a  replica  of  the  transmitted  signal,  and  measure  the  time  delays  between  peaks  of  the  output 
to  determine  relative  positions  of  the  hydrophones.  By  doing  this  with  signals  from  two  (or  more) 
sources  at  different  aspect  angles,  a  two-dimensional  rendering  of  the  array  element  positions  can  be 
obtained. 

The  first  practical  problem  was  obtaining  a  suitable  replica.  There  are  three  possible  options: 

a.  Mathematically  generated  signal, 

b.  Recorded  replica, 

c.  Using  one  of  the  received  channels  as  a  replica. 
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A  mathematically  generated  signal  is  not  a  very  good  representation  of  the  signal  as  received  by 
the  array,  since  by  the  time  it  is  received  it  has  been  distorted  by  the  projection  transducer,  the  trans¬ 
mission  through  a  few  hundred  meters  of  seawater,  and  (minimally)  by  the  receiving  hydrophone. 

A  recording  of  the  test  signals  was  taken  in  a  calibration  pool.  This  high-fidelity  recording  yields 
a  replica  that  is  transformed  by  the  transfer  function  of  the  projector.  However,  it  was  recorded  at  a 
distance  of  only  0.5  meter. 

The  drawback  of  using  one  of  the  received  channels  as  a  replica  is  the  replica  will  then  be  tainted 
by  any  noise  present  at  that  channel  at  the  time  of  the  recording.  Since  the  signal-to-noise  ratio 
(SNR)  is  about  0  dB  in  some  cases,  this  is  a  significant  concern.  Additionally,  appropriate  window¬ 
ing  of  the  signal  is  difficult  in  a  case  of  little  or  no  signal  excess. 

We  investigated  using  each  of  the  three  methods.  The  criteria  for  choosing  the  preferred  method 
is  its  ability  to  produce  a  sharp  and  distinct  peak  in  the  cross-correlation  with  secondary  peaks  mini¬ 
mized  relative  to  the  main  peak.  When  appropriately  windowed  in  time,  the  use  of  one  of  the  chan¬ 
nels  as  a  replica  produced  the  best  results.  The  recorded  signal  and  the  mathematically  pure  signal 
were  both  inferior. 

A  two-step  correlation  was  implemented  to  automate  the  process  of  choosing  a  window  in  time 
that  contains  the  signal  and  to  choose  the  “best”  of  the  received  channels  as  a  replica.  First,  the 
time-reversed  ideal  chirp  (r)  was  convolved  with  each  column  of  B.  The  column  of  B  that  produced 
the  “best  match”  (highest  peak)  was  chosen  as  the  replica  and  then  convolved  with  the  entire  matrix. 
This  “best”  column  of  B  was  windowed  in  time  to  0.3125  second  (1024  samples)  centered  at  the 
location  of  the  highest  peak.  The  result  of  this  cross-correlation  operation  is  called  D. 

Cjojo  =  max(r  *  B) 

b  =  Bj;j_jo_size(r)  to  io-size(r)+1024,jo 

D  =  B*b 

4,  Find  position  of  peaks. 

The  relevant  information  in  the  matrix  D  (output  of  the  replica  correlation)  is  simply  the  location 
of  the  peak  of  each  column.  This  corresponds  to  the  delay  from  the  transmission  of  the  pulse  to  the 
arrival  at  each  hydrophone.  The  units  are  samples,  with  each  sample  representing  305  psec. 

d  =  i:Djj=max(Dj),  Vj 

5,  Interpolate  over  dead  hydrophones. 

Each  of  the  hydrophones  that  were  zeroed  out  in  step  2  gets  a  delay  value  equal  to  the  arithmetic 
mean  of  the  two  adjacent  hydrophones. 

dj=(dj+i  +  dj_i)/2,  j=8, 23, 27, 28, 29, 38. 

6,  Normalize  delays  relative  to  first  hydrophone  in  array. 

The  high-frequency  array  consists  of  hydrophones  3  through  40.  The  uniform  low-frequency 
array  consists  of  hydrophones  44  through  58.  These  will  be  considered  separately,  so  these  two  sub¬ 
sets  of  d  are  separated  out,  and  then  normalized  so  the  first  element  is  0.  Thus,  the  entire  vector 
takes  on  the  meaning  “delay  relative  to  first  hydrophone.” 
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h  =  dj:j  =3  to  40 
I  =  Clj:j=44  to  58 
hj  =  hj-hi,  Vj 
lj  =  lj-ll,  Vj 

At  this  point,  consider  the  anticipated  response  of  the  row  vectors  h  and  1  (high-  and  low-fre¬ 
quency  arrays).  First,  suppose  the  arrays  are  perfectly  straight.  (In  fact,  if  this  were  the  case,  this 
exercise  would  be  unnecessary,  but  it  does  provide  an  adequate  first-order  approximation.)  If  the 
distance  between  the  sound  source  and  the  array  is  large,  then  the  wavefront  impinging  on  the  array 
can  be  approximated  as  a  plane  wave.  This  way,  the  delay  between  any  two  equally  spaced  hydro¬ 
phones  is  the  same.  In  the  case  where  the  source  is  located  along  the  line  of  the  array,  the  delays 
between  arrivals  will  be  maximized,  and  if  the  source  is  broadside  to  the  array,  the  delays  would  be 
identically  zero. 

More  specifically,  if  the  source  is  located  at  endfire  (along  the  line  of)  the  array,  the  predicted 
delays  can  be  derived  as  follows.  The  spacing  of  the  hydrophones  is  equal  to  half  a  wavelength 
(kll). 


m  =  c/2fa  =  cTa/2, 

where  c  is  the  speed  of  propagation  in  the  medium,  generally  assumed  to  be  1500  m/sec,  fa  is  the 
frequency  of  the  array,  and  Ta  is  the  reciprocal  l/fg.  The  period  of  the  array  Ta  in  number  of  samples 
is  simply  the  ratio  of  the  sample  rate  fg  to  the  array  frequency  fa. 

Ta=fs/fa- 

Thus,  the  time  delay  (in  number  of  samples)  between  elements  is 

?i/2c  =  fs/2fa  =  2.56  for  fa  =  640  Hz, 

=  10.24  for  fa  =  160  Hz. 

7,  Find  first  order  polynomial  fit  of  the  delay  data. 

Find  the  straight-line  fit  that  best  approximates  the  delay  data  in  a  least  squares  sense.  If  j  is  the 
hydrophone  (column)  number,  then  we  want  to  find 

Ph  =  [Ph1  Phz]  :  ^j(Ph*[j  1]'  -  hj)^  is  minimized, 
for  the  high-frequency  array,  and  similarly  for  the  low-frequency  array. 

Pi  =  [Pl1  P12]  :  ^j(Pl*[j  1]'  -  Ij)^  is  minimized. 

Note  that  in  the  ideal  case  of  a  source  at  endfire  to  a  perfectly  straight  array  in  the  direction  of 
hydrophone  number  1,  the  polynomials  will  be  Ph  =  [2.56  -2.56],  and  P|  =  [10.24  -10.24]. 

8,  Convert  slope  of  polynomial  to  angle  of  source. 

The  slope  of  the  best-fit  polynomial  can  be  used  to  determine  the  angle  from  which  the  wavefront 
originated.  The  slope  (pi)  is  physically  constrained  to  be  less  than  lfs/2fal.  The  angle  can  be  found 
as: 
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0  =±((fs/2fa-pi)/(fs/fa))7r 

=  ±((2.56  -  Phi)/5.1 2);i,  for  the  high-frequency  array 
=  ±((1 0.24  —  Pii)/20. 48)71,  for  the  low-frequency  array 

Although  the  actual  value  of  the  angle  is  of  interest,  the  first  term  of  the  polynomial,  pi,  is  impor¬ 
tant  for  purposes  of  visualization.  If  each  element  of  the  delay  vector  is  divided  by  pi,  and  the  result 
is  plotted  against  hydrophone  number  on  the  abscissa,  then  the  slope  is  normalized,  and  recordings 
taken  from  different  angles  relative  to  the  array  can  be  overlaid  and  directly  compared. 


COMPARISON  OF  RESULTS 

Table  1  shows  a  summary  of  each  recording  in  the  data  set  that  was  initially  considered.  No 
ground  truth  knowledge  exists  regarding  the  actual  exact  positions  of  the  arrays.  However,  it  is 
known  that  motion  was  negligible,  based  on  the  lack  of  acoustic  interference  that  would  have  accom¬ 
panied  any  motion  of  the  acoustic  elements.  Thus,  the  best  comparison  that  can  be  made  is  one  of 
consistency  of  the  measurements  from  one  trial  to  the  next.  The  reader  is  also  invited  to  compare 
results  from  an  alternate  method  of  element  location  in  Williams  (1994). 

Table  1 .  Summary  of  data  files  used  in  analysis. 


Date 

Name/Time 

CHIRP 

Location 

9/14/94 

G201400 

HF 

BETWEEN 

9/14/94 

G201500 

HF 

OFF  HF ARRAY 

9/14/94 

G201800 

HF 

BETWEEN 

9/14/94 

G201900 

HF 

OFF  HF ARRAY 

9/14/94 

G202100 

HF 

OFF  HF ARRAY 

9/14/94 

G203100 

HF 

OFF  HF ARRAY 

9/14/94 

G204040 

LF 

BETWEEN 

9/14/94 

G204440 

LF 

BETWEEN 

9/14/94 

G204540 

LF 

OFF  HF ARRAY 

9/14/94 

G204840 

LF 

BETWEEN 

9/14/94 

G220100 

HF 

OFF  LF  ARRAY 

9/14/94 

G220140 

LF 

OFF  LF  ARRAY 

9/14/94 

G220500 

HF 

OFF  LF  ARRAY 

9/14/94 

G220540 

LF 

OFF  LF  ARRAY 

9/14/94 

G220900 

HF 

OFF  LF  ARRAY 

Table  2  shows  the  results  of  the  procedure  described  in  the  previous  section  applied  to  the  appli¬ 
cable  data.  Not  all  of  the  chirps  produced  meaningful  measurements  on  both  arrays.  Possibly  some 
of  these  were  missed  due  to  a  timing  problem. 
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Table  2.  Best  fit  polynomial  delays  and  computed  angle  of  sound  sources. 


Name/Time 

Location 

Array 

P1 

ABS(Angle) 

G201400 

BETWEEN 

HE 

-2.4668 

177 

G201400 

BETWEEN 

LF 

9.3643 

8 

G201500 

OFF  HE  ARRAY 

HF 

-0.344 

102 

G201500 

OFF  HE  ARRAY 

LF 

9.1857 

9 

G201800 

BETWEEN 

HF 

-2.1695 

166 

G201800 

BETWEEN 

LF 

9.5 

7 

G201900 

OFF  HE  ARRAY 

HF 

-0.3485 

102 

G201900 

OFF  HE  ARRAY 

LF 

9.3643 

8 

G202100 

OFF  HE  ARRAY 

LF 

9.7929 

4 

G204040 

BETWEEN 

HF 

0.1347 

85 

G204440 

BETWEEN 

LF 

-2.8857 

115 

G220100 

OFF  LF  ARRAY 

HF 

-1.9013 

157 

G220500 

OFF  LF  ARRAY 

LF 

-9.1429 

170 

G220540 

OFF  LF  ARRAY 

LF 

-9.0321 

169 

G220900 

OFF  LF  ARRAY 

LF 

-2.975 

116 

Figures  5  through  10  show  relative  delay  versus  phone  number  for  the  high-frequency  array,  and 
figures  11  through  19  show  relative  delay  versus  phone  number  for  the  low-frequency  array.  The 
title  on  these  plots  refer  to  the  time  the  data  trial  was  recorded.  Refer  back  to  figure  4  for  a  layout  of 
the  arrays  with  the  approximate  source  locations. 

If  the  slope  is  normalized  as  was  described  in  step  8  of  the  procedure  section,  all  trials  on  a  given 
array  can  be  viewed  on  a  common  set  of  axes  for  comparison.  Figure  20  shows  the  six  data  trials  on 
the  high-frequency  array,  and  figure  21  shows  the  nine  data  trials  on  the  low-frequency  array. 

DIRECTION  OF  SOURCE 

The  overwhelming  conclusion  based  on  direction  of  the  source  is  that  a  direction  closer  to  endfire 
(either  end)  gives  better  results.  A  measurement  can  only  be  made  to  one  sample  period,  and  this 
quantization  is  most  detrimental  when  the  arrival  times  are  all  close  together.  G204040  is  the  closest 
to  broadside,  whereas  G201400  is  closest  to  endfire  on  the  high-frequency  array.  Figures  5  and  9 
make  this  point  clear. 

The  last  column  in  table  2  shows  the  calculated  angle  of  arrival  of  the  chirp.  Recall  that  since 
these  are  based  on  a  perfect  linear  array,  they  are  subject  to  a  sign  ambiguity.  They  are  consistent 
with  the  actual  deployed  position  of  the  sound  source  (which  is  not  known  to  any  great  precision)  in 
all  cases  except  one  —  G204040. 

Another  interesting  observation  with  regard  to  G204040  is  shown  in  figure  22.  This  shows  a  sub¬ 
set  of  the  lines  in  figure  20  (G201500,  G201900,  and  G204040).  As  seen  in  table  2,  all  are  nearly 
broadside,  but  the  fact  that  G204040  appears  as  a  mirror  image  of  the  other  two  indicates  that  it  most 
likely  originated  from  a  southwesterly  direction,  while  the  others  came  from  the  northeast.  The  pro¬ 
found  implication  of  this  fact  is  that  the  process  works  well  even  when  applied  to  signals  other  than 
the  intended  ones !  This  will  be  further  addressed  in  a  later  section. 
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Figure  8.  G201 900,  HF  array. 
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Figure  13.  G201800,  LF  array. 


Figure  14.  G201900,  LF  array. 
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Figure  15.  G202100,  LF  array. 


Figure  16.  G204440,  LF  array. 
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Figure  17.  G220500,  LF  array. 


Figure  18.  G220540,  LF  array. 
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Figure  19.  G22900,  LF  array. 


Figure  20.  HF  array,  six  trials. 
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15 


Figure  21 .  LF  array,  nine  trials. 


Figure  22.  HF  array,  G201 500  G201 900  G204040. 


19 


Since  the  normalizing  process  amounts  to  rotating  the  best-fit  line  to  a  slope  of  1 ,  a  comparison 
can  be  made  regarding  how  close  the  data  is  to  that  line.  Figure  23  shows  the  phone-by-phone 
difference  between  the  slope-normalized  data  and  a  line  of  slope  1  for  the  high-frequency  array. 
Figure  24  is  a  corresponding  look  at  the  low-frequency  array. 

Figure  25  is  identical  to  figure  23  except  the  differences  from  trial  G204040  are  inverted  to  reflect 
the  presumption  that  the  arrival  was  from  the  other  side.  This  shows  considerable  agreement  with 
the  remainder  of  the  data  trials  on  this  array.  This  information  is  used  to  “correct”  the  normalized 
slope  depiction  of  G204040  (flipping  the  deviations  to  the  other  side  of  the  best-fit  line),  and  is 
shown  in  figure  26. 

ARRAY  SPACING  RELATIVE  TO  SAMPLE  PERIOD 

A  comparison  of  the  data  on  the  two  arrays,  with  spacing  differing  by  two  octaves,  can  be  used  to 
draw  more  general  conclusions.  The  current  sample  rate,  when  used  with  the  160-Hz  array,  is 
roughly  equivalent  to  using  a  sample  rate  2  octaves  higher  with  the  640-Hz  array.  Clearly  the  quan¬ 
tization  of  the  arrival  times  is  less  of  a  problem  when  fg/fa  is  larger.  Compare,  for  example,  figures  7 
and  17,  which  are  recorded  from  similar  angles. 

More  precisely,  consider  the  consistency  of  the  measurements  for  the  two  arrays.  The  mean  of 
the  slope-normalized  data  for  array  #1  (G204040  flipped)  is  shown  in  figure  27.  The  mean  of  the 
slope-normalized  data  for  array  #2  is  shown  in  figure  28.  As  a  measure  of  consistency,  the  devi¬ 
ations  from  the  mean  for  the  two  arrays  are  plotted  in  figures  29  and  30.  The  standard  deviation  of 
the  measurements  on  array  #1  is  1.545.  On  array  #2,  the  standard  deviation  is  0.372.  Note  that  the 
ratio  of  these  two  values  is  practically  identical  to  the  ratio  of  the  frequencies  of  the  two  arrays. 

BANDWIDTH  OF  PULSE 

The  results  from  the  high-bandwidth  and  low-bandwidth  pulses  were  compared.  Average  vari¬ 
ance  was  0.71  for  the  high-bandwidth  pulse  and  2.36  for  the  low-bandwidth  pulse.  However,  much 
of  the  difference  is  attributed  to  one  extreme  outlier  (G204040)  of  the  three  trials  using  the  low- 
bandwidth  pulse.  Since  the  effect  of  bearing  relative  to  the  source  is  already  demonstrated,  and  the 
few  samples  available  are  not  uniformly  distributed  across  bearing,  no  conclusion  should  be  drawn 
regarding  bandwidth  of  the  pulse. 

EXTRANEOUS  SOUNDS  AS  LOCALIZATION  SOURCES 

As  was  noted  in  the  discussion  above,  trial  G204040  exhibited  a  useful  signal  from  a  direction 
other  than  where  the  known  sources  were  located.  This  was  seen  on  a  few  additional  trials  as  well. 

Figure  3 1  shows  the  output  of  the  replica  correlation  step  of  the  procedure.  This  example  is  for 
G201400,  and  shows  a  good  response  for  all  phones  on  both  arrays  over  a  relatively  short  range  of 
delays.  (The  absolute  numbers  associated  with  the  delay  are  not  meaningful;  the  time  window  is 
selected  to  present  the  relevant  data  at  a  reasonable  level  of  detail.)  This  figure  shows  the  nature  of 
the  response  from  the  intended  signal. 

Contrast  this  figure  with  figure  32,  which  shows  just  the  15  phones  of  array  #2  on  trial  G220500. 
Notice  that  two  signals  are  arriving  nearly  simultaneously  from  different  directions.  If  either  one 
existed  without  the  other,  it  would  apparently  be  sufficient  to  perform  element  localization.  This 
would  seem  to  indicate  that  extraneous  signals  could  be  used  for  the  purpose  of  element  localization. 
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Figure  23.  HF  array,  deviation  from  linear. 
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Figure  24.  LF  array,  deviation  from  linear. 
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Figure  25.  HF  array,  deviation  from  linear  2. 


Figure  26.  G204040,  slopw  normalized,  flipped. 
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Figure  27.  HF  array,  slopw  normalized  mean. 
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Figure  28.  LF  array,  slope  normalized  mean. 
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Figure  30.  LF  array,  deviation  from  mean. 
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Figure  31.  G201400,  Correlation  with  signal  from  channel  6. 


Figure  32.  G220500,  Correlation  with  signal  from  channel  48. 


Figure  33  shows  a  longer  time  window  of  the  same  data  trial  as  figure  32,  as  well  as  displaying  all 
58  channels.  The  high  coherence  between  the  arrays  that  was  exhibited  in  figure  31  (and  others 
within  a  few  minutes  of  it)  is  no  longer  seen.  It  seems  that  the  conditions  deteriorated  in  the  nearly 
2  hours  between  these  trials.  Figure  34  shows  the  correlation  with  a  time-window  of  channel  37  that 
most  likely  includes  the  signal.  Thus,  it  is  seen  that  the  coherence  within  each  array  remains  good, 
although  the  coherence  between  them  is  poor.  Figure  35  looks  very  similar,  but  it  uses  a  different 
and  disjoint  time-window  from  channel  2  as  a  replica;  thus,  it  is  believed  to  be  a  response  due  to  an 
extraneous  signal.  This  suggests  that  the  signal  and  noise  emanate  from  the  same  direction.  One 
scenario  consistent  with  this  fact  is  the  possibility  that  engine  noise  or  cavitation  from  the  boat, 
which  was  used  to  deploy  the  sources,  is  acting  as  the  extraneous  sound  source. 

Figure  36  presents  a  depiction  of  how  good  the  noise  coherence  is  across  array  #1.  This  is  the 
same  data  as  in  figure  35,  but  the  viewing  angle  has  been  altered  to  present  a  side  view  so  the  relative 
heights  of  the  peaks  can  be  seen. 

Another  example  of  two  signals  arriving  from  different  angles  at  nearly  the  same  time  is  seen  in 
figure  37.  The  signal  in  this  case  is  nearly  broadside,  while  the  extraneous  response  is  close  to  end- 
fire.  Using  a  later  time-window  of  the  same  channel  as  a  replica,  the  noise  is  still  seen,  but  the  sig¬ 
nal  no  longer  exists.  See  figure  38. 

In  the  same  pattern  as  figure  29,  the  deviation  of  the  measurments  from  the  unintended  signals 
from  previously  calculated  mean  of  the  controlled  signals  is  plotted  in  figure  39.  The  agreement 
with  the  other  measurements  is  excellent  These  two  trials  exhibit  a  standard  deviation  of  only  0.482. 
Thus,  the  signals  we  acquired  by  chance  seem  to  work  every  bit  as  well  as  the  signals  we  tried  hard 
to  create. 


CONCLUSION 

The  concept  of  processing  employed  was  consistent  enough  to  warrant  further  interest.  The 
method  is  not  reliable  enough  that  a  single  measurement  provides  sufficient  confidence  to  be  used  as 
corrections  to  a  beamformer.  However,  averages  of  about  half  a  dozen  measurements  appear  to  have 
considerable  consistency. 

The  deviations  between  measurements  are  attributed  largely  to  quantization  of  the  signal.  Further 
analysis  could  validate  the  possibility  of  interpolating  the  existing  data  to  a  higher  sample  rate.  Thus, 
the  data  may  not  need  to  be  recorded  at  a  higher  sample  rate,  just  interpolated  prior  to  processing. 

The  conclusion  that  was  drawn  about  bearing  from  source  to  array  (endfire  is  preferred  to  broadside) 
is  also  based  on  time  quantization. 

Extraneous  sources  of  acoustic  energy  were  shown  to  be  of  comparable  utility  in  localization  of 
the  array  elements.  This  also  warrants  further  investigation  when  the  array  is  redeployed.  Con¬ 
trolled  testing  of  broadband  propulsion  noise  at  various  ranges  and  bearings  from  the  array  would 
provide  data  necessary  to  decide  if  such  a  crude  sound  source  is  acceptable. 
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Figure  33.  G220500,  Correlation  with  signal  from  channel  48. 


Figure  34.  G220500,  Correlation  with  signal  from  channel  37. 
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Figure  37.  G2031 00,  Correlation  with  signal  from  channel  6. 
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Figure  38.  G2031 00,  Correlation  with  noise  from  channel  6. 


Figure  39.  HF  array,  deviation  from  mean,  noise  measurements. 

NOTATION 

Matrices  are  denoted  by  bold-faced  capital  letters,  e.g.,  A.  Indices  to  the  rows  and  columns  of  a 
matrix  are  i  and  j  respectively. 

Vectors  (either  row  or  column  vectors)  are  represented  by  bold-faced  lower-case  letters,  e.g.,  b. 
Consistent  with  the  matrix  notation,  elements  of  a  row  vector  are  subscripted  with  j,  and  elements  of 
a  column  vector  are  subscripted  with  i.  If  a  vector  is  a  subset  of  a  matrix,  the  letter  used  will  be  the 
same,  e.g.,  Ci  is  a  row  of  C. 

Scalars  are  plain  text,  e.g.,  f  or  T.  If  a  scalar  is  an  element  of  a  matrix,  it  will  be  upper-case  with 
two  subscripts  (Dy),  and  if  it  is  an  element  of  a  vector  it  will  be  lower-case  and  have  a  single  sub¬ 
script  (fi). 
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